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The paper deals with a next-to-leading order analysis of the radiative B — » A"* 7 decay. Working 
in a PQCD approach, we compute the correction to the essential form factor, coming from a single 
gluon exchange with the spectator. We investigate the supersymmetry effects on the branching ratio 
and direct CP asymmetry and constrain the squark mixing parameter (523) r „■ 



I. INTRODUCTION 

After the Cabibbo-favoured b — > sj mode was first 
reported, in 1993, by CLEO II [1] and updated in 1995 
[2], the exclusive radiative decays, B — > K*^y and B — ► 
fry, as well as the inclusive ones, B — > X s r ( n'y, have 
become main targets for both experimental and theo- 
retical investigations. The exclusive modes, which are 
easier to be experimentally investigated [3, 4, 5], but 
less theoretically clear, have been worked out in dif- 
ferent approaches. For example, the spin symmetry 
for heavy quarks combined with wave function mod- 
els [6, 7] or the heavy quark effective theory when 
both b and s are heavy [8] have been used. Also, per- 
turbative QCD (PQCD) formalisms, introduced for 
exclusive nonleptonic heavy-to-light transitions, have 
been extended to account for the radiative decays. Re- 
cently, detailed analyses of B — > K*^ and B — > ^7, in 
the next-to-leading order (NLO), with the inclusion of 
hard spectator and vertex corrections, have been per- 
formed [9-11] and a consistent treatment, based on 
a new factorization formula, has been proposed [12]. 
Besides an independent determination of the [V^/Vy 
ratio, the b — > sj decays are suitable for studying 
the viability of SUSY extensions of the SM, in view 
of flavour changing neutral currents (FCNC) and CP 
tests, and for imposing constraints on the supersym- 
metric benchmark scenarios [13, 14]. 

The aim of the present paper is to analyse the 
B — > if*7 decay, in the minimal supersymmetric SM 
(MSSM) context. First, at next-to-leading order, 
we compute the hard-spectator correction to the 
essential form factor. In this respect, we employ the 
PQCD approach developed by Szczepaniak et al. for 



decays dominated by tree diagrams [15] and later 
extended to penguin processes [16]. In order to fit the 
Br estimation with data and since large CP asym- 
metries, which are expected to be measured more 
precisely in the near future, at the hadron machines, 
seam to have a new physics origin, we extend our 
analysis beyond the SM. In this respect, we make 
use of the mass insertion method and include, in 
the Wilson coefficients gluino-mediated FCNC 

contributions. Finally, the Br data and the BaBar 
allowed range for direct CP asymmetry are used to 
constrain the squark mixing parameter (S^) LR - 

II. NLO CORRECTION TO FORM FACTOR 

The effective Hamiltonian describing the B — > K*-f 
radiative decay is given by [9, 10] 

H = ^| A p [C 7 7 + dOf + CsOs] , (1) 

where A p = V p bV* s , with p summed over u and c, and 
C\, C7, Cg are the effective Wilson coefficients at /i = 
mj,. The hadronic matrix elements of the four-fermion 
operator and of the electromagnetic and chromagnetic 
penguin operators 

&[ = (^(1-75)^(^(1-75)6) , 

° 8 = [^(1 + 75)^6^, (2) 

possess the general Lorentz decomposition: 



Insert PSN Here 



2 International Conference on Hadron Physics TROIA '07, Canakkale, Turkey, 30 August - 03 September 2007 



(K*\s^b\B) 



2iV{g 2 ) 
m B + m K . 



c- „c* v P a pP 



(K*\s W b\B) = 2m K ,A (q 2 )^-q fl + A 1 (q 2 )(m B +m K ^ 



-A 2 (q 2 ) 



m B + m K * 



(P x P ^ ( m B- m K-) „ ' 



(K*\so-^q"b\B) = 2T 1 (q 2 )e^ a0 e*»P K ,P B <, 

(K*\sc-^ 5 q"b\B) = -tT 2 (q 2 )\(m 2 B -m 2 K ,)e;-(e*q)(P B + P K ,) fi 



-iT 3 (q 2 )(e*q) 



% m% - m 2 K , 



(Pb+Pk*) l 



(3) 



(4) 



where q^ is the the momentum of the photon and e v is 
the K* 4- vector polarization. The form factors are not 
known from first principles and this imprecise knowl- 
edge is a major source of mismatch between theory 
and data as well as between different theoretical esti- 
mations [9]. 

In the heavy quark limit, to;, Aqcd, by neglect- 
ing the corrections of order l/mj, and a s , one has the 
following relation among the form factors [9] 

mB v(o) = mB + mK * A l{ o) 

m B +m K , m B 

= T 1 (0)=T 2 (0) = F K »(0) (5) 

This is broken when one includes QCD radiative cor- 
rections coming from vertex renormalization and hard 
gluon exchange with the spectator. At order a s , the 
form factors encode strong interaction effects by re- 
ceiving an additive correction from hard spectator in- 
teraction [11]. We recommend [9, 10, 12] for detailed 
analyses of both factorizable and nonfactorizable ver- 
tex and hard-spectator contributions, involving the 



operators O7, Os and penguin-type diagrams of 0\. 
However, it has been stated that factorization holds, 
at large recoil and leading order in 1/mj,, and quanti- 
tative tests for proving QCD factorization at the level 
of power corrections have been provided [17]. 

For a consistent treatment of radiative decays, at 
next-to-leading order in QCD, a novel factorization 
formula have been proposed in [12]. In this approach, 
the hadronic matrix elements in (1) are written in 
terms of the essential form factor, which describes the 
long-distance dynamics and is a nonperturbative ob- 
ject, and of the hard-scattering kernels, T- and T/ 7 , 
including the perturbative short-distance interactions, 
as 



(K*7\Oi\B) = \Fk* (0)T/ + 4> B ® T* 1 ® <pK'~\ - V, 

(6) 

where r] is the photon polarization. When the domi- 
nant contribution comes from O7, we use (4) to write 
down the decay amplitude as 



= —7= A. 



G F x em b (n) 



V2 P 2n 2 
and consequently the branching ratio reads 



CV {fj) F k ,(0) \e^ a ^e*"P K * P P B -i(P K , q) {r 1 e* )+ i{e*q) ( V P K , ) 



(7) 



Br LO = r B G >\ X ^< m | (1 _ z 2 f\C 7 (m b )\ 2 \F K , (0)| 2 , 



32tt 4 



(8) 



where z — mK* /m B . At next-to-leading order in a s , 
one has to consider, in (6), the contributions to the 
hard scattering kernels T/ coming from the operators 



0\ and 8 . These have been evaluated in [12] and 
bring (7) to the expression 
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Gf em t ()i) 

A V2 P 2n 2 



47T 



where C F = (N 2 - 1)/(2N), N = 3, and 
Gi(s) = 



(9) 



833 _ 20wr + Stt 2 gV2 + 2 ^ + 3Q _ ^ 2 _ ^3 _ 36C(3) + (3g + ^ _ ^ ^ 



+ (3 + 6i7r) In 2 s + In 3 s] s + - [18 + 2tt 2 - 2i7r 3 + (12 - 6tt 2 ) Ins + 6mt In 2 s + In 3 s] s 2 



1 



G 8 = — 



27 

n _ 2^ 
3 9 



-9 + 112mt - 14tt 2 + (182 - 48m) In s - 126 In 2 s] s 3 , 



2i?r 
~3~ 



(10) 



g 



g 



xP 



yP K 




FIG. 1: The Feynman contributing diagrams in the hard 
scattering amplitude T M . The gluon and photon are re- 
spectively represented by dotted and dashed lines. 



with s c = m 2 /m 2 and fi = rut,. 

Going further, wc add factorizable NLO hard- 
spectator corrections, to the form factor Fk*(0). For 
a single gluon exchanged with the spectator (see Fig. 
1), we extend the PQCD approach, developed by 
Szczepaniak et al. for heavy-to-light transitions dom- 
inated by tree diagrams [15], to the so-called penguin 
processes. 

For the operator O7, one has to evaluate the fol- 
lowing trace over spin, flavor and color indices, and 
integrate over momentum fractions [16] 



Tji = Tr 



+ Tr 



7 /, , x v fa + m b a Ag 2 

(pK'V^i 1 +75)? ,„2 _ ^2 "fa0B7 Jy2 



K~ m b 



,2n 



. (11) 



where Q 2 w —(1 — x)(l — y)m B . The B meson wave 
function 

<pB = ^<PB(x)(PB+mBh5 (12) 



contains a strongly peaked distribution amplitude, 
around a = X B /m B w 0.072, for A B = 0.38. The 
if* is described by the wave function 



fjc- 
12 



<PK'{y)pK*i, 



(13) 



where the light-cone distribution amplitude, ipK*{y), 
has the following expansion in Gegenbauer polynomi- 
als [18] 



<Pk-(v) = 6y(l-y)[l + af* cf /2) (2y-l 
+ afcl'\2y-l) + ...}, 



(14) 



C 3/2 (u) = 3u, 
(m 6 ) = 0.18 ± 0.05, and af (m b ) = 0.03 ± 0.03. 



Ct /2 {u) = (3/2) (5u 2 - 1), 



with 
af 

Performing the calculations in (11) and using the form 
factors decomposition (4), we identify the spectator 
contribution to the essential form factor as 



F sp (a) 



ffs fj3_fK* 

9 nifiAfl 



dy 



(2-y) 
(i-y) 2 



<PK*(y) , 



(15) 

where the K*-m&ss has been neglected. Since we have 
introduced a cut-off in y — > 1, the form factor correc- 
tion (15) depends on the peaking parameter a and this 
is a main uncertainty in our calculations. For the fol- 
lowing input values: a s (p — Q 2 ) w 0.38, fs — 0.180 
GeV, f£. = 0.185 GeV and a = 0.072, we get 



F sp (0.072) = 0.1475 



(16) 



Introducing, in (9), the total form factor Fk*(0) + 
F s p(0.072) = 0.38 + 0.1475, where we have used the 
prediction from the QCD sum rules analysis Fk* (0) = 
0.38 ± 0.06, the branching ratio gets significantly en- 
hanced to the value Br NLO = 6.97 x 10~ 5 . This 
is comparable to the average theoretical prediction 
(7.5 ±0.3) x 10~ 5 [9, 11, 12], but is above the ex- 
perimental data: 
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( (3.83 ± 0.62 ± 0.22) x (BaBar [3]) 

Br(B+ -» K*+j) = I (3.76±g"|| ± 0.28) x lO" 5 (CLEO [4]) 

{ (3.89 ±0.93 ± 0.41) x 10~ 5 (Belle [5]) 

( (4.23 ±0.40 ±0.22) x 10~ 5 (BaBar [3]) 

Br(B° -» K*°j) = I (4.55±°'g ± 0.34) x 10" 5 (CLEO[4]) 

{ (4.96 ± 0.67 ± 0.45) x 10~ 5 (Belle [5]) 



whose world average value, over the and B° decay 
modes, is 

Br exp (B ± -> K*^) = (4.22 ± 0.28) x 10~ 5 (17) 

For the direct CP asymmetry, 

_ T(B -> K*-y) T(B -> K*j) 
aCP ~ T(B — > if*7) + T(B — > K*j) ' (18) 

which is predicted by the SM to be acp < |0.005|, the 
BaBar and CLEO data are: a C p = -0.044 ± 0.076 ± 
0.012 (BaBar [3]) and a C p = 0.08±0.13±0.03 (CLEO 
[4], for the sum of neutral and charged B — > K*j de- 
cays). Even though the data look consistent with the 
SM result within la, it is too early to draw a definitive 



conclusion because of the large errors. With BaBar 
and Belle large data samples, there is hope for pre- 
cise measurements of large CP asymmetries and this 
makes room for new physics. Moreover, it has been 
stated that CP asymmetries larger than few percent 
would have a dominantly supersymmetric origin [19]. 



III. BRANCHING RATIO AND DIRECT CP 
ASYMMETRY WITHIN MSSM 

Following this idea, let us analyse the B — ► K*-f 
decay in the MSSM context. Using the mass insertion 
approximation, [20], we incorporate, in the Wilson co- 
efficients C 7 and C 8 , the FCNC SUSY contributions 



where 



n SUSY , M \ 

W ( m susy) = 
C s (Msusy) = 



V2n 



G F (V ub V u * s + V cb V*)mj 



m b 



F (x) 
G (x) 



G F {V ub V: s + V cb V* s )m 2 - v Ii)LR m b 



4a; 



9(1 -a;) 4 
x 



[l + Ax - 5x 2 + Ax ln(x) + 2x 2 \n{x)] 



3(1 -x) 4 



[22 - 20a; - 2a; 2 + 16a; ln(x) - x 2 ln(x) + 9 ln(a 



(19) 



(20) 



In (20), x = m\ ]/m~ is expressed in terms of the gluino Wilson coefficients, encoding the new physics, become 



mass, m,g, and an average squark mass, iriq. We un- 
derline that, in the expressions of Cy § SY (Msusy), we 
have kept only the left-right squark mixing parame- 
ter (£33) LR = (Afcg) /m? since, being proportional to 
the large factor mg/m b , it has a significant numerical 
impact on the branching ratio value. The quantities 
Ab s are the off-diagonal terms in the sfermion mass 
matrices, connecting the flavours b and s along the 
sfermion propagators. In these assumptions, the total 



C* ota '[x,«5] = C 7 (m b ) + C 7 SUSY (m b ) , 
Cl otal [x,5] = C & (m b ) + Ci USY (m b ), (21) 

where C 7 g SY (m b ) have been evolved from Msusy = 
m g down to the /j, = m b scale, using the relations [21] 

Ci USY (m b ) = V Ci USY (m s ), 
C 7 SUSY (m b ) = V 2 C 7 SUSY (m s ) 

+ l(v-V 2 )Ci USY (m- g ), (22) 



Insert PSN Here 



International Conference on Hadron Physics TROIA '07, Canakkale, Turkey, 30 August - 03 September 2007 5 



with 

n = (a s (mj)/« s (m t )) 2/21 (a s (m t )/a s (m i ,)) 2/23 

(23) 

Finally, putting everything together, we replace, in 
(9), the Wilson coefficients CV and C 8 respectively 



by C\°\ f [x, 5} and Cl° tal [x,5], the form factor F K " (0) 
by .F K *(0) + F sp (a) and, consequently, the branching 
ratio is 



Br total = Br SM+SUSY = Tg 

x A 



G\arn\ 3 



32tt 4 
Cp 



m^il-z 2 ) 6 \F K * (0) + F s P(a) 
2 



C\° tal \x, S] + (CiG? + Cl otal [x, 5}G 8 ) 



(24) 



For a given x and 5 = pe lip , the total branching ratio 
(24) is depending on three free parameters: a, p, p. 
One can notice that, by including the hard gluon con- 
tributions as a correction to Fk*(0), the direct CP 
asymmetry parameter, (18), is free of the uncertain- 
ties coming from the form factors. 

In what it concerns the gluino, as its pair production 
cross section has large cancellations in the e + e~ an- 
nihilation, there is hope that the laser-backscattcring 
photons will provide a precise gluino mass determina- 
tion [22] . For a wide range of squark masses, a gluino 
mass of 540 GeV may be measured, with a precision 
of at least ±2 ... 5, at the multi-TeV linear collider at 
CERN. 

In the next coming discussion, we use the following 
input parameters: mb(mb) — 4.2 GeV, a — 1/137, 
\V tb V t *\ = 0.0396 ± 0.002, t b o = (1.546 ± 0.018) ps, 
uiq = 500 GeV and x has the values xi — 0.3, xo = 
(540/500) 2 and x g = 3 (where 1(g) comes from m g 
less (greater) than m,q). 

In figure 2, we represent the contour plots on which 
the Br total , (24), is equal to the world average data 
(17) and the BaBar constraint [3] 

- 0.17 < a CP < 0.082 (25) 

When {p, if} e [0, 0.03] x [— tt/2, tt/2], we get for (17) 
three dashed lines, with increasing thickness, as x goes 
from xi to x g . Correspondingly, for acp, we get three 
pairs of solid curves: the lower ones, for acp — —0.17, 
and the upper ones, for acp — 0.082. These solid con- 
tours close inside the values of direct CP asymmetry 
which do not agree with (25). One is able to con- 
strain (S23) LR to the segments of the -Br-plots out- 
side the solid contours, for each x. For example, when 
m g > niq, all negative phases, with suitable p's, can 
accommodate both (25) and (17). Also, p is con- 
strained by the experimental data on the branching 
ratio of B -> X s -/ to p < 0.016 [20, 23]. 

In figure 3, we represent the Br total (in units 
of 10~ 5 ) and acp (in units of 0.1), with respec- 
tively dashed and solid lines, as functions of ip, for 



9- 





A N 


N / / J 


i 







0.005 0.01 0.015 0.02 0.025 
P 



FIG. 2: Contour plots of total branching ratio, in units 
of 1CT 5 , fitting the world average data, (17), (the dashed 
lines) and the BaBar constraint (25) (the solid lines), as 
functions of p and ip. The thickness of contours is increas- 
ing as x is taking the values: x = {0.3, 1.16, 3}. The solid 
curves close inside the values of acp which disagree with 
the constraint (25). 



x = xq. As p takes the following values: p e 
{0.005, 0.01, 0.012}, we get three pairs of curves, with 
increasing thickness. The horizontal dashed band 
corresponds to the data 3.76 x 10~ 5 < Br exp < 
4.68 x 10~ 5 (in 90% C.L.), while the horizontal solid 
band stands for the constraint (25). We notice that 
one should avoid the region ip e [— 7r/16, 7r/16] on the 
p w 0.01 plot, since the CP asymmetry parameter 
drops quickly, from positive to negative values much 
outside the constraint (24). 

Finally, let us perform a numerical analysis, for 
x = .To, and increasing p, starting with p — 0.005. As 
p € [-87T/16, -4tt/16] U [37I-/16, 7tt/16], the Br total 
and the direct CP asymmetry are inside the ranges 
10 5 x Br total e [8.3, 3.3] and [3.1, 8.1], accommo- 
dating data and other theoretical models predic- 
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FIG. 3: The total branching ratio, in units of 10 -5 , (the 
dashed lines) and 10 x acp (the solid lines), as functions 
of if, for x — xq. The thickness of plots increases as p 
is: p = {0.005, 0.01, 0.012}. The horizontal dashed band 
corresponds to the world average branching ratio with 90% 
C.L. and the horizontal solid one is for the constraint (25). 



tions, and, respectively, acp £ [—0.054, —0.093] U 
[0.096, 0.062]. When p goes to bigger values, the 
two ip ranges, constrained by the allowed branch- 
ing ratios, get closer and acp moves toward much 
bigger values. For example, for p — 0.01 and 
ip £ [-671-/I6, -471-/16] U [37r/16, 5tt/16], one gets 
10 5 x Br total £ [8.1, 3.6] and [3.2, 7.62] and, respec- 
tively, a CP £ [-0.1, -0.16] U [0.21, 0.12]. We notice 
that only the negative <p- values lead to acp inside the 
BaBar constraint (25). For p close to the upper limit, 
p = 0.015, and ip £ [-4tt/16, -2tt/16]U [tt/16, 3tt/16], 



the values 10 5 x Br total £ [7.8, 3.6] and [3.3, 7.3] are 
compatible with a measurable acp ~ ±0.12. 



IV. CONCLUSIONS 

In the present paper, we have analysed the ra- 
diative B — > K*j decay in the MSSM framework. 
First, we have used the PQCD approach, developed 
by Szczepaniak et al. [15] and extended to "penguin" 
processes [16], to compute the hard-spectator contri- 
bution, F sp (a), to the essential form factor Fk» (0). 
For the peaking parameter in the B wave function 
a = 0.072 and Fk* (0) = 0.38, the branching ratio be- 
comes Br NLO = 6.97 x 10~ 5 , which is above the ex- 
perimental data. In order to reach an agreement with 
data and find large CP asymmetries, which hopefully 
will be measured with a better precision in the near 
future, we extend our analyses by including, in the 
Wilson coefficients CV^, the SUSY contributions com- 
ing from squark mixing parameter (#23) LR — P e%lf '■ 
Consequently, the total branching ratio depends, be- 
sides a, on three (SUSY) parameters: x, p, ip, while 
acp is free of the uncertainties coming from the cor- 
rected form factor. Using the graphs displayed in 
figures 2 and 3, one is able to find ranges for the 
mass insertion parameter (^23)^- As an example, 
for x = (540/500) 2 , the world average branching ra- 
tio, Br exp = 4.22 x 10~ 5 , can be accommodated for 
{p, ip} = {0.005, -|f} or {0.01, --[§}. The cor- 
responding asymmetries, acp — —0.085 and respec- 
tively acp = —0.147, are inside the BaBar constraint 
(25). 
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